Phase transition properties of a finite ferroelectric superlattice 

from the transverse Ising model 



Xiao-Guang Wang 1 ' 2 * Ning-Ning Liu 2 Shao-Hua Pan 1,2 and Guo-Zhen Yang 1,2 
1. China Center of Advanced Science and Technology (World Laboratory), 
P.O.Box 8730, Beijing 100080, People's Republic of China 
2. Laboratory of Optical Physics, Institute of Physics, 
Chinese Academy of Sciences, Beijing 100080, People's Republic of China 

(February 1, 2008) 

Abstract 

We consider a finite ferroelectric superlattice in which the elementary unit cell 
is made up of I atomic layers of type A and n atomic layers of type B. Based 
on the transverse Ising model we examine the phase transition properties of 
the ferroelectric superlattice. Using the transfer matrix method we derive 
the equation for the Curie temperature of the superlattice. Numerical results 
are given for the dependence of the Curie temperature on the thickness and 
exchange constants of the superlattice. 
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I. INTRODUCTION 



Possibly because of the great difficulty of growing well characterized samples, experimen- 
tal studies of ferroelectric superlattices have been published only in recent years(Iijima et al. 
1992; Tsurumi et al. 1994; Wiener-Avnear 1994; Tabata, Tanaka and Kawai 1994; Kanno 
et al. 1996; Zhao T. et al. 1999). Some exploratory theoretical work on ferroelectric super- 
lattices has appeared (Tilley 1988; Schwenk, Fishman and Schwabl 1988; Schwenk, Fishman 
and Schwabl 1990). Their starting point is the Ginzburg-Laudau phenomenological theory. 

On the microscopic level, the transverse Ising model (TIM)(de Gennes 1963; Binder 
1987; Tilley and Zeks 1984; Cottam, Tilley and Zeks 1984)was used to study infinite ferro- 
electric superlattices under mean field theory (Qu, Zhong and Zhang 1994; Qu, Zhong and 
Zhang 1995; Zhong and Smith 1998 ) or effective field theory(Zhou and Yang). From the 
experimental point of view the TIM is a valuable model because of its possible applications, 
for example in studies of hydrogen bonded ferroelectrics(de Gennes 1963), cooperative Jahn- 
Teller systems(Elliot et al. 1971) and strongly anisotropic magnetic materials in a transverse 
field(Wang and Cooper 1968). The reviews of Blinc and Zeks(1972) and Stinchcombe(1973) 
give more details about possible applications of the TIM. 

In the present paper, we consider a finite ferroelectric superlattice in which the elementary 
unit cell is made up of I atomic layers of type A and n atomic layers of type B. The mean- 
field approximation is employed and the equation for the Curie temperature is obtained 
by use of the transfer matrix method. We study two models of the superlattice which 
alternate as ABAB...AB (Model I) or ABABA.. .BA (Model II). Numerical results are given 
for the dependence of the Curie temperature on the thickness and exchange constants of the 
superlattice. 
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II. THE CURIE TEMPERATURE 



We start with the TIM(de Gennes 1963; Sy 1993, Qu, Zhong and Zhang 1994; Qu, Zhong 
and Zhang 1995; Zhong and Smith 1998; Bouziane et al. 1999) 

H = — - ^ ^ J ijS z ir S z - r , — ^ QiSjf r , (1) 

(r,r') ir 

where Sf r ,Sf r are the x and z components of the pseudo-spin, are plane indices and 

(r, r') are different sites of the planes, J^- denote the exchange constants. We assume that 
the transverse field fli is dependent only on layer index and consider the interaction between 
neighboring sites. For simplicity, we take Q the same in the superlattice because the main 
qualitative important features result from the difference of J^. 
The spin average (Si) , obtained from the mean field theory 

where Hi(Q,0,J2j Jij(Sj)) is the mean field acting on the ith spin, ks is the Boltzman 
constant and T is the temperature. 

At a temperature close and below the Curie temperature, (Sf) and (Sf) are small, \Hi\ rs 
Q, equation (2) can be approximated as 

{s ' } = ^ tallh (^r)^ J «< s ^ + 2 ( J -'+'< s J + i> + J -.--i< s ti))] W 

Here zq and z are the numbers of nearest neighbors in a certain plane and between successive 
planes respectively. 

Let us rewrite Eq.(4) in matrix form in analogy with the reference(Barnas 1992) 
lm i+1 \ M { m ,\ 

with Mi as the transfer matrix defined by 

((r — zqJh)/ (zJii+i) —Ji,i-i/ Ji,i+i \ 
o )■ (6) 

where rrii = (Sf) and r = 2Q/(zJj i j + i) coth[Q/ (2 ksT)}. 



We consider a ferroelectric superlattice which alternates as ABAB...AB . In each ele- 
mentary unit AB, there are / atomic layers of type A and n atomic layers of type B. The 
intralayer exchange constants are given by J A and J B whereas the exchange constants be- 
tween different layers is described by J AB . We assume there are N elementary units and the 
layer index is from to N(l + n) — 1. In this case, the transfer matrix M; reduces to two 
types: 

(X A -1\ (X B -1\ 

Mi oj-Mi «)• (7) 

where X A = t-j'a,^b = r - j B , 3a = z J A /(zJ AB ),j B = z J B / (z J AB ), and t = 
2n/(zJ AB ) coth[Q/(2A; B T)]. 
From Eq.(5), we get 

\m N{l+n) _ 2 ) \m J 

where 

R = M B ~ 1 ( M l A M B ) ^ _ 1 M l A 1 (9) 

is the total transfer matrix. 

From the above equation and the following equations 

m 1 = X A m ,m N{ i +n) _ 2 = X B m N{ i +n) _ u (10) 

we obtain the equation for the Curie temperature of the superlattice as 

RnX A X B + R U X B — R 2 ±X A — R 2 2 = 0. (11) 

Next we consider Model II, the superlattice which alternates as ABA...BA and assume 
that the lattice has N(l + n) + I layers. The total transfer matrix 

S = M^MbR (12) 

and the equation for the Curie temperature is obtained as 

S U X 2 A + (S 12 - S 21 )X A - S 22 = 0. (13) 



For an unimodular matrix M, the n-th power of M can be linearized as(Yariv 1992; 
Wang, Pan and Yang 1999) 

M n = U n M — U n -il, (14) 

where I is the unit matrix, U n = (A™ — A")/(A + — A_), and A± are the two eigenvalues of 
the matrix M. 

Using Eq.(14), we obtain 



M\ = Ei Ma - E^J, (15) 
M n B = F n M B - F n _J, (16) 

where E x = (a l + - a l _)/(a+ - a_), F n = (#» -02 )/(/?+- P-), a± = (X A ± yj X\ - 4)/2 and 
/3± = (X B ± ^X 2 B -A)/2. Then from Eqs.(15) and (16), the matrix M^M™ in Eq.(9) can be 
written explicitly as 



M AB = M\M n B = 

(EiX a — Ei-i) (F n X b — E n -i) — EiF n 

Ei (F n X b — E n -i) — E X -\F n 

The trace of the matrix M AB is 

tr = {EiX a - £,_!) (F„X 6 - K-i) - 2E t F n + E^E^. (18) 

Since det (Mab) = 1, the eigenvalues of the matrix Ma B is 7± = (tr ± \Jtr 2 — 4)/2. Then 
using Eq.(14), we get 

M^ 1 = Gat-iM as - GW_ 2 J, (19) 

where = (7^ - 7- )/(7+ - 7-)- 

Using Eq.(15)-(19), we can express the total transfer matrix R and S in terms of 
Xa, X b , EiF n , and G^. We can get an explicit expression for equations (11) and (13) for the 



— (E[X a — Ei-i) F n + E x E n -\ 
—EiF n + Ei-\E n -\ 



(17) 
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Curie temperature by substituting the matrix elements of R and S into Eq.(ll) and (13), 
respectively. The results are tedious, we only give numerical results below. 

Fig.l gives the dependence of the reduced Curie temperature tc against the reduced 
exchange constant ]a in model I and II. The Curie temperature increases with increase of 
2a- It is clear that the Curie temperatures in model II are larger than those in model I. The 
reason is that the superlattice in model II is thicker than that in model I. The fact that the 
Curie temperature increases with the increase of j A can also be seen in Fig.2. Fig.2 shows the 
dependence of the reduced Curie temperature tc against the reduced exchange constant 2 a 
for different ui in model I. The transverse field causes a reduction of the Curie temperature. 
In other words, the Curie temperature decreases with increase of u>. 

Fig. 3 shows the dependence of the Curie temperature on the number of elementary units 
N in model I. to m the figure is the Curie temperature of the corresponding infinite superlat- 
tice. The Curie temperature of the infinite superlattice can be determined from the following 
equation(Wang, Pan and Yang 1999) 

trace(M^M^) = 2. (20) 

The Curie temperature of a finite superlattice is always less than that of a corresponding 
infinite superlattice, and it increases with the increase of the number of elementary units N 
to approach asymptotically to t for large values of N. 

III. CONCLUSION 

In conclusion, we have studied the phase transition properties of a finite ferroelectric 
superlattice in which the elementary unit cell is made up of / atomic layers of type A and 
n atomic layers of type B. By the transfer matrix method we derived the equation for the 
Curie temperature of the superlattice. Numerical results are given for the dependence of the 
Curie temperature on the thickness and exchange constants of the superlattice. The method 
proposed here can be applied to the finite superlattice in which each elementary unit cell is 
made up of many types of materials and the atomic layers of each type can be arbitrary. The 
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finite superlattice is more realistic than the infinite superlattice in experiments. We hope 
that the present work will have relevance to some future experiments. 



Captions: 

Fig.l, The dependence of the reduced Curie temperature tc against the reduced exchange 
constant ja in model I and II. The parameters js — 1, I — n — N — 2, and u = 0.5. 

Fig. 2, The dependence of the reduced Curie temperature tc against the reduced exchange 
constant ja for different uj in model I. The parameters js = l,and / = n = N = 2. 

Fig. 3, The dependence of the Curie temperature on the number of elementary units iV 
in model I. The parameters ja = 1-2, Jb — 1,1 — n — 2, and u> = 0.5. 
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